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We present an application of the Lorentz model in which fits to vibrational spectra or a Kramers 
Kronig analysis are employed along with several useful formalisms to quantify microscopic charge 
in unoriented (powdered) materials. The conditions under which these techniques can be employed 
are discussed, and we analyze the vibrational response of a layered transition metal dichalcogenide 
and its nanoscale analog to illustrate the utility of this approach. 
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I. INTRODUCTION 

Born and local (or ionic) effective charge are well- 
known quantities with which to assess chemical bonding 
and polarization in a material. ^'^ A large local effective 
charge indicates a highly ionic system (for instance, 0.8 in 
NaCl), a medium-sized value points toward an intermedi- 
ate bonding case (0.4 in GaAs), and a small value is asso- 
ciated with a covalent material (0 in Si).''''* On the other 
hand, Born effective charge describes the static and dy- 
namic polarizations.^ From the optical properties point of 
view, large longitudinal optic-transverse optic (LO-TO) 
splittings are well-known characteristics of polarizable 
compounds and associated with substantial Born effec- 
tive charges. This is because LO-TO splitting is directly 
proportional to charge within Born's original formalismii 
Using the LO-TO splitting and high frequency dielec- 
tric constant, it is relatively straightforward to quan- 
tify chemical bonding from optical measurements of high 
quality single crystal samples and compare the extracted 
value(s) with first principles calculationsj^i^ Unfortu- 
nately, there are many instances when single crystals of 
a bulk material are unavailable, either because they can 
not be grown or are not of sufficient size or quality for 
optical measurements. At the same time, nanoparticles, 
nanotubes, and alloys (or composite mixtures) present 
scientifically compelling problems^^ i-'^^'^^i^^'^'^i"'^'*'^^ where 
optical measurements on powdered materials are the only 
option, a drawback that complicates the situation but 
does not diminish the desirability of obtaining quantita- 
tive Born and local effective charge data. 

Following Born and Szigetijii^ we present an applica- 
tion of the Lorentz model in which fits to vibrational 
spectra or a Kramers Kronig analysis are employed along 
with several useful formalisms to quantify microscopic 
charge in unoriented (powdered) materials, assuming 
that the effects of ionic displacement and atomic polar- 
izability are superimposable. We demonstrate that this 
technique can be used to assess chemical bonding and 
local strain under certain conditions, a development that 
advances the field of nanoscience and, at the same time, 
retains many attractive features of optical spectroscopy 



and the traditional Lorentz model. This paper is or- 
ganized as follows. Sections 11 and 111 review the situ- 
ation for single and multiple coUinear oscillators. Sec- 
tions IV and V present how this approach must be mod- 
ified for a randomly oriented (powder) sample. Section 
VI illustrates the use of this technique to quantify Born 
and local (or ionic) charge in powdered 2H-M0S2 and 
its nanoscale analog. Comparison of our analysis of the 
powdered 2H-M0S2 data with that of the single crystal 
shows that this approach accounts almost perfectly for 
sample orientationJ^iilii^ The extension to assess size 
effects in nanoparticle samples demonstrates its utilityJ^ 
Our objective is to clearly present the framework, use- 
ful equations, and an application of this method. The 
Appendix provides alternate frameworks of the model. 



II. COLLINEAR OSCILLATOR 

A. A Review of the General Formalism 

Starting with the Lorentz model for bound charge^ii^ 
one has: 



mx + m-fx + mujQX = qE, 



(1) 



where x is the relative displacement of positive and nega- 
tive ions, q is the effective ionic charge, m is the reduced 
mass, E is the electric field, 7 is the damping parameter 
and luq is the oscillator frequency. This model is ap- 
plicable to vibrational modes of isolated oscillators such 
as those in gas phase molecules. Due to dipole-dipole 
interactions between oscillators in a solid, so called de- 
polarization effects must be considered. Here, E has to 
be replaced by the microscopic field Emic 



P 

E ^ E„u „ = E + 7]—, 



(2) 



where E^ic is the microscopic field that is really acting on 
the oscillator, E is the macroscopic field that should be 
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used to calculate the dielectric constant, P is the polar- 
ization, and rj is the depolarization factor which depends 
on the topological arrangement of the oscillators. With 
this substitution, Eq. ([1]) becomes: 



mx + rwyx + muj^x = q{E + 77 — ). 

Co 



(3) 



The overall goal is to find the dielectric constant which 
relates E and P. However, in Eq. ([3]), there are 3 un- 
known variables: x, E and P. One more equation is 
needed to find the relationship between E and P. This 
additional information will come from the definition of 
P. Assuming the polarization P is a linear combination 
of ionic contributions Pi = which comes from rela- 
tive displacements of ions, and electronic contributions 
Pg = ^^Emic which is due to the distortion of electron 
clouds, one has 



P = P^ 

xq 



Pe 



V 



V 



where a and V are the polarizability and volume of the 
oscillator, respectively. Substituting the expression for 
E^ic in Eq- ©I one has the second important equation: 



P 



xq 

v 



V 



E 



P 



or 



l-'qa/V^V V ' 



(4) 



Taken together, Eqs. ([3]) and (j4]) are enough for dealing 
with single crystal problems (where there is translational 
symmetry). Next we discuss the dielectric response. 

To obtain the expression for the dielectric constant 
which relates P and E^ we must eliminate P From Eq. 
©. Plugging Eq. (g]) into Eq. Q, one has 



mx + mjx + itluiqX = q[E + 



7] 



,xq epa 



eo l-'qa/V^V 



mx + m-fx + [mtjQ 



q^T] 



V 



qE 



epV{\-'qa/Vy l-'qa/V 
For X = xoe"'"', the solution is 



m{l — ria/V) 



E 



- meoV{l-nc./V) \ " ^ " 



(5) 



For simplicity, we define: 



,2 _ 



^ meQV{l~T]a/V)' 
Using Eq. we can write down the polarization: 
1 ,xq Cna . 

1 r q m{l-7ja/V)^ 



eoa/V 



l-r]a/V 



E 



E 



Therefore, the dielectric constant is 



(6) 



e = 1 



1 + 



P 

a/V 



1 



r]a/V 

„2 



(7) 



We can immediately see that the Lorentz model in a 
solid is modified compared with that of an isolated os- 
cillator due to the depolarization effect and the polar- 
izability (when a = and ry = 0, Eq. ([7|) reduces to 
the case of an isolated oscillator. ^^). We end this sec- 
tion by summarizing several useful definitions and ex- 
pressions that connect measurable quantities (left-hand 
side) to microscopic parameters (right-hand side). These 
include the high-frequency dielectric constant, oscillator 
strength, and TO phonon frequency. 



e(oo) 
A 



a/V 



1 - -qa/V 

„2 



eomV{l - rjOilVf 



2 2 2 



q^rj 



eQmV{l — rjOi/V) 



(8b) 



(8c) 



Clearly, one can extract e(oo), A, and uj^o from the ex- 
perimental vibrational spectra using an oscillator fit or 
Kramers Kronig analysis and use these quantities to ex- 
tract the microscopic parametersj-^ In other words, once 
e(oo) and A are measured, one can derive a and q using 
Eqs. (Hil) and (l8b|. Note that 



S . 



■"TO 



is the dimensionless oscillator strength. 



B. Evaluating the Effective Charges 

Ionic effective charge (also called local effective charge) 
q is an important quantity because it quantifies the ion- 
icity. It is distinct from the Born effective charge as dis- 
cussed below. The most straightforward way to extract 
ionic effective charge from oscillator strength is to use 
Eq. (US): 



= eomVil - i^a/VfA. 



(9) 



Thus, quantitative information about bond cova- 
lency/ionicity can be extracted from a knowledge of os- 
cillator strength and polarizability. 

In the literature, we often find Born effective charge 
defined asi^i 



1% 



eomV 

which, as shown in the Appendix, is equivalent to: 

„2 



9% 
eomV 



A 



eomV{l - rja/V)^' 



(10) 



(11) 



Then the dielectric constant and other definitions can 
be expanded as: 



e = 1 



j_ 



eamjV{l ~ Vjaj/Vf (wg - wf_^-) - ^2 - i^^ 



= l + ^[£,(oo)-l]+^A, 



,(14) 



where 



£j(oo) = 1 



^2 



eomjVil-r]jaj/V)' 



(15a) 
(15b) 

(15c) 



Therefore Born and local ionic charge are related as 



1% = 



{l~r]a/V)^ 



(12) 



Note that Born effective charge takes into account the 
combined contributions of ionic displacement, electron 
cloud deformation, and depolarization effects, whereas 
ionic effective charge only represents the charge of the 
ions. 



Note that polarizability represents the high frequency 
dielectric response of the electron cloud. Therefore, a 
should be labeled according to the polarization direction, 
although for simplicity, we label this quantity with the 
mode index. 



IV. TILTED OSCILLATOR 



III. MULTIPLE COLLINEAR OSCILLATORS 

When there is more than one oscillator (as in most real 
materials), one has to go back to the polarizability (Eq. 
^) and add a mode index j. This yields: 

P = "^/^ E 
' 1 - v,a,/V 

I g| E 

mjV{l — r]jaj/Vy oJq — ^1 j — — i^juo 

(13) 



If the electric field is not perfectly aligned with the 
direction of motion of a certain mode, the observed os- 
cillator strength will be reduced from its intrinsic value. 
This can be easily understood by considering the extreme 
case: when the light is polarized perpendicular to a par- 
ticular mode, this mode will not contribute to the os- 
cillator strength at all. Thus, if one directly employs 
the formulas for local ionic and Born effective charge in 
Section Hi Bl without modification, the results will be un- 
physical. This is because a tilted oscillator provides only 
a component of the total intrinsic oscillator strength. In- 
stead, we must go back to Section III Al and rederive a set 
of formulas that take orientation into account. 

We can employ a modified version of Eq. ([5]) to account 
for the effect of a tilted oscillator: 



and 
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Hence, 

p{e) 



.xq eoa 



l-ria/V'V V 
eoa/V 



Ecos{e)] 



l-Tja/V 



Ecos{e) 



Ecos{0) 



and 



ei9) = 1 + 
= 1 + 



P{9) cos{e) 

a/Vcos^{6) 
1 - ria/V 

„2 



0082(6*) 



eomV{l — rja/V)"^ lu^ — ujf — lu^ — ijuj 



(18) 



Here, 6 is the angle between the electric field and the 
direction in which oscillator intensity is maximum. If a 
measurement is done on a sample with a distribution of 
orientations {{0}), the result is appropriately averaged 
as: 

e - {e{e)} 

^ a/V{cos^{e)) 



1-ria/V 

„2 



(COS2(0)} 



eomV{l — rja/V)^ uJq — loI — uj^ — i^uj 



(19) 



In this case, the observed parameters are related to the 
microscopic parameters as: 



e{oo) = 1 + 



a/V{cos'^{e)) 
l-ria/V 



A 



e{)mV{l -~ rja/V)' 



(20a) 



■(cos2(6')) (20b) 



(20c) 



° eomVil -rja/V)' 
Here, the brackets () indicate directional averaging. 

V. MULTIPLE TILTED OSCILLATORS 

Most isotropic samples of real materials have several 
vibrational modes. We can write down the dielectric con- 
stant for the case of multiple tilted oscillators by combin- 
ing Eqs. in]) and P7)) : 



(21) 



eomjV{l - TljUjlVy (w^ - Luf -) - uj"^ ~ i-fjOj' 



(22) 



and 



eim) = l+5^[e,(oo)-l]cos2(0,)+^A, 



cos 



Note that the {6j} are related to an oscillator's polariza- 
tion direction. Therefore, the number of independent Oj 
may be less than the number of modes. Hence, 



e = {e{{e,})) 

= l + ^[£,(^)-l](cos2(0,)) 



e(oo) + YAj 



(COS2(0,)) 



(23) 



•"TO.j 



where the observed (apparent) parameters are: 



eioo) = 1 + E^ 



A, 



~2 



J 

l + ^[e,(^)-l](cos2(0,.)) (24a) 
j 

(cos2(0,)) (24b) 



2 2 
, ,2 



eQmjVil-r^jaj/V) 



(24c) 



If all vibrational features can be resolved in frequency 
space, e(oo) and Aj can be determined from a model 
oscillator fit or a Kramers-Kronig analysis. At the 
same time, the rjj can be estimated from the crystal 
structure. 2" However, one can not use the -I- 1 equa- 
tions given by Eqs. (|24all24bl) to solve for 3iV unknowns. 
The latter include a,-, g'j, and {cos'^{6j)). Even if in some 
cases, we know {cos {9j)) (perhaps from an independent 
x-ray measurement), there are still 2A'^ unknowns to de- 
termine from only N + 1 equations. Additional informa- 
tion is needed to constrain the system. 



Despite this limitation, vibrational spectroscopy of un- 
oriented powdered samples can still be an important tool 
for extracting microscopic charge and bonding informa- 
tion. There are two important cases: 

Case 1: The system is simple enough to have + 1 = 
2N . For example, in rocksalts, the crystal symmetry is 
cubic ((cos^(6'j)) is always 1/3), and there is only one 
vibrational mode (iV = 1). Thus powder spectroscopy is 
sufficient to determine all of the microscopic parameters 
for systems such as NaCl or MnO. 

Case 2: Occasionally, some variables are already 
known, say, from another method, sample, or similar 
compound, so that the total number of known variables 
can be reduced to be equal to or less than + 1. Recent 
work on M0S2 nanoparticles provides a good exampleJ^ 
Here, the interplane oscillator strength is identical to that 
of the single crystal. Therefore, it is reasonable to assume 
that the corresponding polarizability and charge are both 
the same for the nanoparticles as they are in the single 
crystal, a coincidence that reduces the number of un- 
knowns and makes the extraction of intra-plane charge 
bonding information possible. We elaborate on the case 
of M0S2 below. 



VI. THE DYNAMICS OF A MODEL 
TRANSITION METAL DICHALCOGENIDE: 
TESTING OUR APPROACH 

In order to test the workability of this approach, we 
elected to investigate a model transition metal dichalco- 
genide. The bulk material, 2H-M0S2, belongs to the 
PQs/mmc space group (Fig. [Tfa)).^^ One consequence of 
this layered architecture is the low-dimensional electronic 
structure which consists of strong bonding within layers 
and weak van der Waals interactions between layersi^^ 
Each M0S2 slab contains a layer of metal centers, sand- 
wiched between two chalcogen layers, with each metal 
atom bonded to six chalcogen atoms in a trigonal pris- 
matic arrangement. There are two infrared active Eiu 
and ^211 vibrational modes f22, Schematic views of these 
displacement patterns are shown in Fig. [TJa). The 
Eiu and symmetries characterize intralayer and in- 
terlayer motions, respectively. We begin with demon- 
strating the self-consistency of the theory by analyz- 
ing the oscillator orientation and predicting the ob- 
served optical parameters. We then extend our technique 
to the chemically-identical but morphologically different 
nanoparticles to illustrate the consequences of finite size, 
strain, and curvature. 

Figure [Ijb) displays a close-up view of the reflectance 
of 2H-MoS2.^"' As expected for a pressed powder sam- 
ple, both the Eiu and modes are clearly observed. 
At normal incidence, the dielectric function is related to 
reflectance as: 



IR Active Modes (a) 
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2u ; 
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FIG. 1: (Color online) (a) Crystal structure of 2H-M0S2,— 
schematic view of the displacement patterns of the infrared- 
active modes, photo of a typical pressed pellet sample, and 
diagram for tilted oscillators, (b) Close-up view of the 300 K 
reflectance spectra of 2H-M0S2 powder. The red solid curve is 
experimental data and the blue dashed line is a flt according 
to Eqs. (|25|l and (|28|) . The insets show the dielectric constant 
and a schematic view of the 2H-M0S2 platelets in a pressed 
pellet sample, (c) Close-up view of the 300 K reflectance of 
IF-M0S2. An identical color scheme and line type is employed. 
The insets show the dielectric response and a high resolution 
TEM image of an IF-M0S2 nanoparticle. 



R{uj) 



(25) 



Equation (P5|) is formally valid for single-bounce re- 
flectance at the interface of 2 semi-infinite media. For 
the case of real materials with finite thickness, sample 
thickness must be sufficient to assume that there is no 
back reflectance. Large attenuation due to a strong res- 
onance is helpful here. 

A fit to the 2II-M0S2 powder data using Eqs. 
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TABLE I: Apparent parameters extracted from an oscillator 
fit to the measured reflectance spectrum of 2H-M0S2 powder, 
and the intrinsic parameters from a single crystal samplcfi^ 







?(oo) " 








0.fl4 


10.3 


0.20 


15.2 


A2u 


0.0036 


0.030 


6.2 



"Apparent parameters obtained from an oscillator fitting analysis 
of the measured powder spectrum. 

''Intrinsic parameters obtained from a fit to the measured re- 
flectance of a single crystal sampleiH 



and ((^ (Fig. [TJb)) allows us to extract £(00) = 10.3, 
Se^^ =0.114 and Sa^^ = 0.0036 (Table I). We refer to 
these values as "apparent oscillator parameters" because 
they are obtained from direct fits to the measured powder 
spectrum. They are distinct from the "intrinsic param- 
eters" that are useful for evaluation of Born and ionic 
charge. 

The apparent and intrinsic oscillator parameters are 
related according to Eqs. ()24at - I24c|) by the oscillator 
orientation. For M0S2, one has: 

+ [eA2„(oo)-l](cos2(0^,J}z; (26a) 
Se,^ = SE,Acos\eE,J)v (26b) 
Sa,^ = SA,Acos''{eA,J)v. (26c) 

Here, v — 0.7 is the relative density of the unori- 
ented pressed pellet compared with that of the single 
crysta L^*^'^^ A benchmark is now needed to related the 
apparent and intrinsic oscillator parameters. 

Fortunately, the polarized infrared reflectance of a 2H- 
M0S2 single crystal has been studied by Wieting et al}"^ 
Fits to the reflectance yield e(oo)£;j^ — 15.2, e(oo)^2„ — 
6.2, Se,^ = 0.20, and Sa^^ = 0.03 (Table |Tj). These 
intrinsic parameters are appropriate benchmarks for the 
our pellet samples, because they are made of ^.m sized 
powders, for which the surface effect can be ignored as a 
good approximation. 

In Eqs. (|26al - [^5c)) . the only unknowns are {cos^ {9 EiJ)) 
and (cos^(6'yi2„)}- Fig. [Hb) shows a schematic view of the 
pressed pellet of 2H- powder. Since the polarizations of 
the two modes are orthogonal, one has Oe^^ + 0yi2u = §• 
Hence, {cofi^ [9 e^^)) + (cos^(0a2u)) — 1; meaning that 
there is only one unknown, say, {cos^ {9 e^^)) ■ Using Eqs. 
P6bll26c|) . we find two independent values as 0.81 and 
0.83 in good agreement with each other. To further check 
the self-consistency, we calculated e'(oo) using Eq. (|26ap 
using an average value of {cos'^{9ei^)) = 0.82, yielding 
e^(cx)) = 9.8, in excellent agreement with that obtained 
by direct fitting techniques {^{oo) — 10.3). In addi- 
tion to confirming the validity of our approach, this self- 
consistency also shows that the pressed powder sample 
of 2H-M0S2 will have the same Born effective charge as 
the single crystal ; — which, of course, it must. The 



dielectric response ei was calculated using the intrinsic 
parameters of Ref. [la and is plotted in the inset of Fig. 
[TJb). The dispersive response is typical of an anisotropic 
insulator with two vibrational modes. 

The availability of ^30 - 70 nm average diameter 
nested M0S2 nanoparticles provides an opportunity to 
investigate the impact of finite length scale effects on 
chemical bonding. Figure [T]^c) displays a close-up view 
of the measured far infrared reflectance along with an 
oscillator fit. The apparent parameters obtained from 
this fit can be scaled toward a set of intrinsic oscillator 
parameters using the orientation and density corrections 
outlined above^i^ From an analysis of these intrinsic os- 
cillator strengths and high frequency dielectric constants 
along with mode frequencies, we can extract Born and 
local effective charges for the nanoparticles.^^ In the in- 
tralayer direction, we find that the Born effective charge 
of the nanoparticles is 0.69 e in the intralayer direction, 
significantly lower than that of the layered bulk (1.11 e). 
Here, e is the charge of an electron. We attribute this 
difference to structural strain (and resulting change in in- 
tralayer polarizability) in the nanoparticles The Born 
effective charge of the nanoparticles remains unchanged 
in the interlayer direction (0.52 e). The dielectric con- 
stant was again calculated using intrinsic parameters'^ 
and is plotted in the inset of Fig. [^c). Clearly, the 
dispersive response of the nanoparticles is different from 
that of the bulk in the intralayer direction. They arc the 
same in the interlayer direction. 

Extension of Born and local (or ionic) charge con- 
cepts to nanomatcrials is an important advance because 
most are not (and may never be) available in an oriented 
form.^iii^iiiii^iii^iii Indecd, emerging mechanical and tri- 
bological applications of nanoscale M0S2 require bulk 
quantities of powder with careful size-shape control but 
no orientational control. At the same time, the relation- 
ship between engineering properties such as solid state 
lubrication^^'^^'^'* and the microscopic aspects of charge 
and bonding is an open and interesting question that de- 
serves further study. 



VII. CONCLUSION 

We present an application of the Lorentz model in 
which fits to vibrational spectra or a Kramers Kronig 
analysis of the refiectance are employed along with sev- 
eral useful formalisms to quantify microscopic charge and 
polarizability in unoriented (powdered) materials. This 
paper provides a systematic development of the opera- 
tive equations and a discussion of the conditions under 
which such techniques can be employed. We demonstrate 
the workability of our approach by analyzing the vibra- 
tional response of a layered transition metal dichalco- 
genide, and we include an evaluation of Born and local 
(or ionic charge) of its nanoscale analog to illustrate the 
modern utility. The extension to assess size effects ad- 
vances the field of nanoscience and, at the same time. 
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retains many attractive features of optical spectroscopy 
and the traditional Lorentz model. 



Ionic effective charge 
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g2 = ei^mV{l-ria/VfA 

= eiimV{l-ria/Vfujlo[e{Q)~-e{oo)] 



= eomV 



^To[£(0) -g(oo)] 



772 [£(oo) + l/7?-l]2 

• Born effective charge 



(29) 



APPENDIX: OTHER FRAMEWORKS 



(30) 



For applications, it can be convenient to use other 
forms of Eq. ([7]) that are written in terms of parameters 
that are more related to the experimental observations. 
Before stepping into that, we take a careful look at Eq. 
([7]) and note that there are only three independent pa- 
rameters: a, g, and luq. Other three- variable sets offer 
equivalent expressions. Some of these are detailed below. 



1. e(0), e(oo) and luto framework 

This is a very useful framework because e(0), £(oo), 
and ujxo straightforwardly extracted from the 

optical spectra. 



a. Collinear Oscillators 

Equation (|8ap is already very close to employing this 
new set of parameters. Let's define: 



l-rja/V' ' tQmV{l - rja/Vy u!? 



TO 



e{oo) -f 



(27) 



Then, 



A 



mV{l - r]a/Vf 



with parameters e(0), £(c»), and wf^o, one can rewrite 
Eq. ©: 



By definition, Wio, e = for a longitudinal mode. 
Using Eq. (|^ (and ignoring i"fuj), one has 



/ N 9 £(0) — e(oo) 



■"TO 



which gives the Lyddane-Sachs- Teller relation^ 



LO 



■'TO 



m_ 

£(oo) 



(31) 



(32) 



One can write out uj\o terms of microscopic param- 
eters a, q and loq using the Lyddane-Sachs- Teller relation 
as follows: 



-'LO 



^TO 



e{oo) 



,2 ^ ' '^TO 



^TO- 



^TO 



s{oo) 
A 



r(oo) 



q^r] 



eomV{l — rja/V) 

-,2 



q- 



2 , 1 

° eomV 1 + (1 - r^)a/V 



(33) 



Note that t^\o l^'i'gc than Wq, while lo^q smaller than 



e = £(oo) -I- ujIo 2 2 ■ 



(28) 



2. Multiple Collinear Oscillators 



It is a bit tedious, but not so difficult to show that: 



2_ 2 £(oo) + 1/77-1 



£j-(0) -£j-(oo) 

(34) 
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3. Tilted Oscillators 



• Born effective charge 



e = l + [e(oo)-l)](cos"(0)) 



4- Multiple Tilted Oscillators 



(35) 



q% = eQmVe{oo){LulQ - J^^) 



2. Multiple Collmear Oscillators 



(39) 



= l + 5][£,(oo)-l]cos2(^,) 



3 



ej(0) -ej(oo) 



■ cos 



TO J- - ^ - «7j ^^ 



(40) 



5. Tilted Oscillators 



'={e{{ej})) 

= l + ^[e,(oo)-l](cos2(0,)) 



e.j{Q) ~ £j{oo) 



(cos^(0,)) (36) 



1 + [e(oo) - l](cos2(6')) 



+ e(oo) 2 9 



'l^^{cos'{e)) (41) 



B. e(oo), tJro and Luto framework 
1. Collinear oscillators 

Anotfier choice is to use e[oo), lo^q^ ^lo- Us- 
ing the Lyddane-Sachs- Teller relation, Eq. ([5^ . one can 
eliminate e(0), which yields 



e = e(oo) + e(cx))-^ = — 

Ionic effective charge 



(37) 



4- Multiple Tilted Oscillators 



+y:^.-M ,r^"-^:r';:!:, cos^(^.-) (42) 



j 



e^{e{{e,})) 
= l + ^[£,(oo)-l](cos2(0,)) 



eQmV{l~r]a/VfA 

<iomV{l - ■qa/VYe{(X)){uj'lQ - uj'^q) 



r]'^[e{oo) + 1/77 - If 



(38) 



/ 2 2 

£j -2 "^-2 (COS^ {O3)) 



(43) 
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